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Abstract 

We extend a dimensional upper bound on how much an optimal trans- 
port map can degenerate for the quadratic transportation cost, originally 
due to Caffarelli, to cost functions that satisfy the curvature condition of 
Ma, Trudinger, and Wang. 



1 Introduction 

The present paper addresses how much an optimal transport map can degener- 
ate in terms of the dimension of certain sets. In optimal transport theory, one 
considers two probability measures p and p, on domains il and fl, respectively. 
There is a given transportation cost function c : O x O 4 I. One is inter- 
ested in the optimal mapping T, a measurable map that minimizes the total 
transportation cost 

/ c(x, F{x))dp{x) 
Jn 

among all F : CI — >• O with F#p = p. For certain type of cost functions, 
including the distance squared cost on Ricmannian manifolds, the existence 
and uniqueness of T, is well established (see [TJ[l)J[T13[2Hl[iH] ), an d it is well 
known that the map T is associated with a potential function, say it, called 
the Brenier solution, and the map T can be found from the (almost everywhere 
defined) gradient of u. In particular, if c(x, x) = —x-x on fi x f2 G M 71 x M. n , then 
we have Vu(x) for almost every x, where u is a convex function. This function 
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u solves the classical Monge- Ampere equation (at least in an appropriate weak 
sense) , 



det(D 2 u(x)) 



p(x) 



p{Wu{x)) ' 



Under appropriate geometric conditions, especially the Ma-Trudinger-Wang 
condition of the cost function and c-convexity of the domains, and when the 
source and target measures p and p are bounded above and below, the optimal 
map T is continuous, and u is differentiable: see [2lgl[8|fT ^ [2lH25 l [28 l l52 l l33]. 
These geometric conditions are in fact necessary for regularity [25l[28]. Also, 
without both bounds on p and p, regularity can no longer be guaranteed: see 
[36]. 

In this paper, our main interest is in the case when the map T is not contin- 
uous, or equivalently when the potential function u is not differentiable. We use 
the affine dimension of the subdifferential du, to measure the degeneracy of the 
map T . Roughly speaking, this dimension measures how many directions the 
function u will be non-diffcrentiable, equivalently, how much the optimal map 
T can split a point and spread its image in the target domain, and thus how 
strong the discontinuity of T can be. See Definitions 12.41 and 12.51 in Section [2] 

For the Euclidean transportation cost c(x, x) = —x ■ x on K™, it is known due 
to Caffarelli [5] that if the source measure p is bounded from above on f2 and 
the target measure p is bounded from below on its support spt then for each 
point x where du{x) (~1 (spt z/) mt ^ 0, the affine dimension of du{x) is strictly less 
than ^ . A Pogorelov type counterexample, also presented in [5] , shows that this 
bound is sharp. What Caffarelli actually showed is the same dimensional bound 
for the contact set, i.e. the set of points having the same image by the optimal 
map, which is equivalent to the bound on the subdifferential by reversing the 
role of the target and the source. 

We extend the dimension estimate of Caffarelli to the case where the cost 
function satisfies the Ma-Trudinger-Wang (|MTW[) condition [28][32]. The main 
theorem we prove in this paper is the following. The relevant definitions and 
conditions concerning c-convex geometry are given in Section [2] below. The 
method presented here is a geometric approach, and it also provides a new 
proof for the case of the Euclidean quadratic cost function. 

Theorem 1.1. Consider two open sets Q C M and Q C M in Riemannian 
manifolds M and M , and fix two probability measures p on f2, and v on Q 
with bounded supports spt p and spt v. Assume p is absolutely continuous with 
respect to Lebesque measure. Suppose that c satisfies $M\i - $M\i and (jMTWj) 
and that and fi are c-convex with respect to each other. Further assume that 
f2 is strongly c-convex with respect to spt v. Finally, let u be the Brenier solution 
of the optimal transportation problem with cost c from p to v , which satisfies 



|3 c uCE)nspti/|£ < ME\c 



(1) 



for any measurable E C ft and some constant A > 0. 
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Then, for any x G fl such that d c u{x) D (spt^) mt ^ 

ft 

affdim (du(x)) < —. 

where affdim is the affine dimension of a convex set. 

Note that the condition d c u(x) (~l (spt^) mt 7^ in the above theorem is 
necessary, since otherwise we have an easy exampie (see for example, [7]) in the 
Euclidean case c(x,x) = —x ■ x, where spt^ = {x £ R 2 \ \x\ < 1}, sptz^ = {x E 
M 2 I 1 < \x\ < 2}, and the Brenier solution u satisfies du(0) = {x \ \x\ < 1}, 
thus, affdim(9w(0)) = 2. 

We also comment here on some variants of the condition (|MTW|) . The first 



is (MTW+ I, which is for example, satisfied by distance squared cost functions on 
the round sphere [26] and its perturbation and quotients pi HUl[TlHT71ll91l27U34j . 
Another variation is (|NNCC|) which is for example, satisfied by the distance 
squared cost on the Euclidean space and the products of round spheres [TU] (see 
also [12]). Both QMTW+P and (|NNCC|) are strictly stronger versions of (jMTWp . 



and one does not imply the other. See Remark 12.21 for the relevant definitions. 
It is worth mentioning that under (|MTW + P or (jNNCCj) . the result in The- 



orem 11.11 is fairly easy to prove. If one assumes ( MTW+ ) , the conditions in 
Theorem 11.11 above are enough to apply the methods in [25] to obtain that the 
c-subdiffcrcntial is a singleton at each point in the domain. On the other hand, 
assuming (INNCCI) , by using the results of [TTJ Section 6] one can directly extend 
the method used in the proof of the main theorem in [5] to obtain the desired 
estimate. 

However, under (|MTW[) alone, it is not clear how to extend Caffarelli's proof. 
The main difficulty lies in the fact that Caffarelli's proof relies on combining 
the arithmetic-geometric mean inequality with the divergence theorem. In the 
case of a more general cost function c, there is a corresponding positive definite 
matrix that the arithmetic-geometric mean can be applied to, but there seems 
to be no easy way to use the divergence theorem to continue the proof. 

The present paper is organized as follows: In Section [2] we give various 
definitions and some standard conditions related to c-convcx geometry and the 
optimal transportation problem. In Section [3] we first provide a proof of the 
main result, but with the Euclidean cost function c{x, x) = —x-x fTheorcm l3.ip . 
This allows us to highlight the geometric framework of the proof without the 
technical difficulties posed by a nonlinear cost function. In Section 2] we provide 
the preliminary setup for the proof of Theorem 11.11 for a general cost function 
c, and in the following Section [S] we give the actual proof. In Section [X] in 
the Appendix, we demonstrate an example of a cost function (which originally 
appeared in |32j and was communicated to the authors by Neil Trudinger) that 



satisfies (|MTWp . but does not satisfy (|MTW + P nor (jNNCCp 
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2 c- convex geometry 



We suppose that f2 and f2 are two bounded, open domains contained in two 
Riemannian manifolds M and M respectively (the bar does not represent the 
closure of a set, we will use the superscript E cl to denote this instead). We will 
also use x and x to denote points in Q and Cl respectively, while D and D will 
either denote the differential of a function with respect to the i or i variables 
respectively. 

Remark 2.1. Suppose we have fixed a coordinate basis on an n-dimensional 
vector space V. Then, for a point, say p £ V , we will use superscripts to 
denote the coordinate entries, and the notation p' = (p 1 , . . . ,p k ) and p" = 
(p k+1 , ■ . ■ ,p n ) to denote the first k and last n — k coordinate entries with respect 
to this basis respectively. Moreover, by an abuse of notation, we will sometimes 
write x' and x" to denote the n-dimensional vectors (i',0) and (0,x"). 

Also, we will write \v\ and v ■ w to denote the standard Euclidean norm 
and inner product for vectors v,w £ R". When the entries are vectors from 
arbitrary vector spaces (the two vectors may be contained in different spaces), 
it is understood that there are fixed coordinate systems (which will be clear 
from context) on the vector spaces, and we are identifying each vector with the 
point in 1™ that gives its coordinate representation. 

Finally, if V is a vector space, the notation (/, v) will denote the action of 
/ 6 V* on the vector v E V. 

We now state a number of conditions, standard in the theory of optimal 
transportation (see [55] and j32], also [TTHTHIHSJUT] f° r more details on these 
conditions). Much of the notation that follows is similar in exposition to what 
is used in [21] by the second author. 
Regularity of cost: 

ceC 4 (!] cl xO d ). (AO) 

Twist : 

c satisfies condition (|A1[) if the mappings 

x i v Dc(x§ 7 x) . 
x i-> —Dc(x, Xq) 

are injective for each Xq € 51 and for each So G Here, D (resp. D) denotes 
the differential in the x (resp. x) variable. 

We denote the inverses of the above mappings by exp% {-) and exp^ (•). 
Notice that for the cost c(x,x) = —x ■ x on R™ x R™, these mappings are both 
just the identity map. 
Nondegeneracy: 

c satisfies condition (|A2[) if, for any x € f2 and x € O, the linear mapping 
given by 

-DDc(x, x) : T s Cl T_ Dc(x ^ (T*Q) = T*n (A2) 
is invertible (hence, so is its adjoint mapping, —DDc(x, x) : T^Cl — > T x fl). 
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The next condition, originally introduced by Ma, Trudinger and Wang [251 
132] , is crucial in the regularity theory of optimal maps and is actually a necessary 
condition for regularity as shown by Loeper in |25j . 
MTW : 

c satisfies condition (|MTW[) if, for all x G fi, x G fi, and any V G T x £l and 
rj 6T*fl such that (77, V) = 0, we have 

- {cij, m ~ '•.,.- '•' •''•..,.,:•<•''•'• ' /f / (.r. x)V l V 3 n km > 0. (MTW) 

Here all derivatives are with respect to fixed coordinate systems on M and M, 
with regular indices denoting coordinate derivatives of c with respect to the 
x variable, and indices with a bar above denoting coordinate derivatives with 
respect to the x derivative. Also, a pair of raised indices denotes the matrix 
inverse. 

Remark 2.2. There are two common, stronger versions of condition (jMTWI) . 

We say that c satisfies condition ( MTW+ ) if there exists some 5q > such 
that for all x G fi, x G f2, and any V G T X Q and 77 G T*Q such that (77, V) = 0, 
we have 

- {c ijm - c^c f ' s c sm )c p ' k c^(x, x)V l V^ m > <5 |T/| 2 M 2 . (MTW+) 

We say that c satisfies condition (|NNCC[) (see [19]) if for all x G £1, x G Cl, 
and any V G T X Q, and rj G T*Q, we have 

- {c ljm - Cl ^c f > s c sm )cV- k c^\x, x)V l V^ kVl > (NNCC) 

(note that (|NNCC|) removes the "orthogonality condition" of (77, V) =0). 

The condition (|MTW[I looks complicated but, as we will see below, it leads 
to some elegant geometric implications: see Lemma 12.81 and Corolloraries 12.91 

andenni 

We now give the definition of a c-convex function, along with its c-subdiffercntial. 

Definition 2.3. We say that a real valued function u defined on fi is c-convex 
if for each point xo G f2, there exists at least one xq G f2 and Aq € 1 such that 

-c(x ,x ) + A = u(x ), 
— c(x, x ) + A < u(x) 

for all x G £1. Any function of the form — c(-,xq) + \q satisfying the above 
relations is said to be a c- function that is supporting to u from below at xq. 

Definition 2.4. For a scmiconvex function u, its subdifferental at xo is defined 
by 

du{x ) = T*n I u(exp XQ («)) > u(xq) + (v,p) + o(\v\)}, 
where exp a , o is the usual Riemannian exponential map at xq. 
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Definition 2.5. If u is a c-convex function and xq G f2, we define its c- 
subdifferential at xq by 

d c u(xo) : — {x G 5l | — c(-, x) + A is a c-function 

supporting to u from below at Xo, for some A G K}. 

Note there is the immediate inclusion, d c u(x) C exp c x (du(x)) for each If 
E C f2, we write 

<9 C «(-E0 := (J <9 c it(x). 

Remark 2.6. Given any x G f2 and x G f2, we denote -E C Q and E C Ct 
represented in the cotangent spaces at x and x respectively, by 

[£]- : = -Dc(E,x), 
[E] x : = -Dc(x,E). 

Definition 2.7. If E C f2 and i G fi, we say -E is (strongly) c-convex with 
respect to x if [E] s is a (strongly) convex subset of T?tt. 

H E C £1 and E 1 C Cl, we say is (strongly) c-convex with respect to E if E 
is (strongly) c-convex with respect to every x G E. 

Analogous definitions hold with the roles of x and x, and E and E reversed. 

Now, suppose we have fixed a point xo G fi. Then, by (|A1|) and (|A2|) . the 
map p i— > exp% a (p) is a diffcomorphism between and H. We then define 

the modified cost function 

c(p,x) := c{exp% o {p),x) - c(exp c So (p), x ) (2) 

for p G [0]^ and x G 51. Also, given a c-convex function u on f2, we modify it 
by defining 

u(p) := u(exp% a (p)) + c(exp5 (p),x ), (3) 

it is easy to see that u is c-convex, and moreover 

x G d c u{exp% a (p)) x G dcu(p). 

When dealing with these modified functions (and also (jTj?)) ) defined on cotangent 
spaces, we will continue to use the notation D and D to denote the differential 
with respect to the p and p variables respectively. It will be clear from the 
context whether the differentiation is with respect to variables on the original 
domains, or on the cotangent spaces. We note here that with this convention, 

-Dc{p, x) = [-DDc{exp% (p), x )] 1 {-Dc(exp% (p), x) + Dc{exp% Q (p), x )) 

(4) 
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These modified functions will be useful when we carry out the proof of the main 
theorem. 

A key result detailing certain geometric properties of c-convex functions was 
discovered by Loeper [35] for domains in R" , further developed in [T71ITO1I371I5T] , 
and extended to domains in manifolds under suitable conditions. 

Lemma 2.8 (Loeper's maximum principle [25J). Suppose c satisfies (|A0[) - (|A2|) 
and (|MTW[) . Also assume fl is c-convex with respect to ft. Fix xq G £1 and let 
p : [0, 1] — > \fi\ x be a line segment. Then, for any x G fl 

- c(x, exp c X0 (p(t))) + c(x , exp c xa {p{t))) 

< max{-c(x, exp c xo (p(0))) + c(x Q , exp c Xa (p(0))), 

- c(x 7 exp L Xo (p(l))) + c(xo,exp xo (p(l)))}. 

An analogous inequality holds with the roles of and Q reversed. 
This lemma has several important consequences: 

Corollary 2.9. H3 Theorem 3.1] Suppose that c satisfies (JJ0])H|M) and (|MTW|) . 

u is a c-convex function, and fl, Q are c-convex with respect to each other. Then, 
for any xq G ft, the set d c u(xo) is c-convex with respect to xq. 

Corollary 2.10. J2B[ Proposition 44] Suppose that c satisfies (|ATJ|) -([A^ ]) 
and (|MTW[) . u is a c-convex function, and ft is c-convex with respect to Cl. 
Then, for any xq G & and Xq G R, if the function u(-) + c(-, xq) — Ao has a local 
minimum at a point xq G fl, then 

x G d c u(xo). 

The following quite useful consequence of Lemma 12.81 is first observed and 
utilized in [TT1[T3] and [2"2] . 

Lemma 2.11. Suppose that c satisfies (|A0|) (| A2|) and ([MTWP . u is a c-convex 
function, and Cl, Cl are c-convex with respect to each other. Then, for any xq G fl 
and Xq G R, the c-sublevel set 

{x G O u(x) < —c(x, Xq) + Xq} 

is c-convex with respect to xq. In terms of the modified function u (see ([5)1 J. the 
function p i-> u(p) is sublevel set convex in p. 

Definition 2.12. Suppose c satisfies conditions (jA0[) " (|A2[) . and [i and v are 
probability measures defined on f2 and f2 respectively. A c-convex function u is 
a Brenier solution for the optimal transportation problem with cost c from fi to 
v if 

T # n = v, 

where the mapping T{x) is given by 

T(x) := exp c x (Du(x)) for a.e. x G Vt. 
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Remark 2.13. It is well known (e.g. see [2HH3H]) that if H is absolutely con- 
tinuous with respect to volume measure on M, and c satisfies (|A0[) - (|A2[) . there 
exists a unique Brcnicr solution (defined a.e.) up to translation, of the optimal 
transportation problem with cost c from /i to v. 

3 Euclidean case 

In this section, we give a proof for the case c(x, x) = —x ■ x on R™ x R n . (One 
should note that this case is equivalent to c(x,x) = \\x — x\ 2 .) This detour 
serves two purposes. First, our method gives an alternative proof of Caffarelli's 
original result [5]. Second, we can better focus on illustrating key geometric 
ideas in our method in the Euclidean case, since in this special case we can 
bypass many of the obstacles that a nonlinear cost creates. 

First note that for this cost function, c-convexity of a set coincides with the 
usual notion of convexity, and a c-convex function is just a convex function. 
The c-subdifferential d c u(xo) at each point xq € f2, coincides with the ordinary 
subdifferential du(xo) as long as fi contains du(xo). This is the case if, for 
example, C spt;/ and u is Brenier solution from /j to f. If E C R", we will 
write 

du{E) := (J du(x). 

xeE 

Moreover, by taking the Legendre transform twice, u can be extended to a 
convex function on all of R ra with 

du(K") C conv(spt/i), 

where conv(-E) denotes the convex hull of a set E. We denote this extension 
also by u. 

In the following we give an alternative proof of Caffarelli's result [S]: 

Theorem 3.1 (see [5]). Consider probability measures fi and v on l n with 
bounded supports spt \x and spt v. Assume \i is absolutely continuous with respect 
to Lebesque measure. Let u be a convex potential on R™, that is a Brenier 
solution for the optimal transportation problem with cost c(x, x) = —x ■ x from 
to v, which also satisfies 

\du(E) nspti/| £ < A\E\ C (5) 

for any measurable E C R n and some constant A > 0. Then, for any x G R" 
with du{x) n (spt^) int ^ 0, 

Ti 

affdim (du(x)) < — , 

where affdim E is the affine dimension of a convex set E. 
Proof. 



Step 



Suppose that the theorem does not hold. Let 

k := affdim (du(x )) = max {affdim (du(x)) x £ R" and du{x) n (spt i/) int ^ 0} 

for some xo £ M. n with fc > - . The idea of the following proof is to first find a 
point, say x e , and two interrelated geometric shapes Qd in K™ (the domain of 
du) and Qd in R™ (the target of du), depending on a small parameter d. Then 
by using the inequality ([5]) , we draw a contradiction as the parameter d goes to 
zero. 

Step 1 

First, let us find a suitable point x e . Note that by ([5]) we must have k < n — 1. 
By subtracting an affine function, we may assume that the origin is contained 
in both the relative interior of du(xo) and (spt^) lnt , and that u(xo) = (in 
particular, this implies u > 0), and by a rotation we may assume that the affine 
hull of du(xo) is spanned by {ei} k =1 , the standard orthonormal Euclidean basis. 
Now define the (n — /c)-dimcnsional affine set 

A n ~ k := p| {x£R n \(x-x o )-x = 0} = {x £ R n \ x' = 0} (6) 

(see Rcmark l2.1l for the notation x' and x"). Also define the contact set, 

S := {x £ R n | u(x) = 0}, 

which is convex since u > 0. Note that So, which is obviously closed, is in fact 
compact. This is from the boundedness of spt /Lt, the fact that du(x)n(spt v) lnt ^ 
0, and ([5]). For example, this can be shown by applying the boundary- not-to- 
interior lemma [l3l Theorem 5.1 (b)], to a sufficiently large ball (a strongly 
convex set) containing spt fx. Since So cannot intersect the boundary of this 
ball, by connectedness it must remain bounded. 

Now, by boundedness of So, we may find a point x e in dSo with the exterior 
sphere property. By a translation of the domain we may assume x e to be the 
origin, i.e. x e = 0. Since du(xo) contains the fc-dimensional ball {x £ K n | 

< r , x" = 0} for some < r < 1, we also find that 

= u(x) > ±(x - a: ) • -^rr 
\x\ 

for any x £ So and x £ du(xo), from which it follows that 

So c A"~ k . (7) 

Moreover, from the fact that u(0) = 0, and the definitions of A n ~ k and the 
subdifferential, it follows 

du(x ) C du(0), 

and by the maximality of k, we find that du(0) also has affine dimension k. 
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Step 2 

Next, let us construct the geometric shapes Qd and Qd- Recall that is a point 
in 8Sq with the exterior sphere property. Hence, by a rotation in x" , we may 
assume for some Rq > that 5*0 is contained in an (n — fc)-dimensional ball with 
radius Rq, i.e. 

S Q c{xe R" | x' = 0, \x" + R a e n \ < R a }. (8) 

We now define some auxiliary sets. For each < p < it (p will be chosen 
later to be close to 7r/2), define the n-dimensional cone 

JC P := {x G R" | |x|cosp < x n } 

(here, p is the angle of the opening of the cone). Then, for d > and < p < ir, 
define 

W d , P := {x e R" | \x"\ = d} n fCp n Bi(0) d . 

We next determine the dependancc of p on the parameter d, as we wish to 
obtain Wd, P D S*o = 0, or cquivalently, > for all a; £ Wd, p . To this end, 
fix x G Wd,p- Note that if x' ^ 0, we have x £ So by ((HJ. On the other hand if 
i' = we can calculate 

|.t" + i? e„| 2 = \x"\ 2 +R% + 2R x n 

> Rl + d 2 +2R \x\ cos p 
= Rl + d 2 + 2R dcosp, 

thus if we choose cos p > — we obtain \x" + i?oe„| > Rq, hence again x ^ Sq 
by ((SJ. Since for d > sufficiently small, 

fit d \ f d \ d 
cos — I = — sin > , 

we will have > for all a; G Wd. Pd , i.e. Wd lj0(J n So = 0, for the choice 

^ := 2 + 4i^- (9) 
Now, we construct and Qd- First, define 

f(d) := inf u(x) 

VV d,p d 

for each d > sufficiently small. Notice that by the previous paragraph, u > 
on Wrf.p d which is also compact by definition, hence we see that for any d > 
small, 

f(d) > 0. 

Observe that since spt v is bounded and du(M. n ) C conv(spt/i), u is Lipschitz 
and therefore f(d)/d < 1. Then, define the set 

Q d := {x G K" | la/I < /(d), |a/'| < d}n/C Pti . 
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We also define 



IC d :={xeR n | |i"|cos(— ) <x n }, 



a cylinder whose base is an (n — fc)-dimensional cone. Notice that the angle jj^ 
here is complimentary to the angle p d . Using this we define the set 

W d := iC d n{xeR n \ \x'\ < \x"\ = r -^} (10) 
and then Q d as the cone over W d with vertex 0, that is, 

Q d := (J [0,5] 
sew d 

where [0,5:] denotes the line segment (in R n ) connecting to x. 
We can see that this set has the volume bound 



n— k 

-k-1 



(P-*- 1 ~ \Q d \ £ (11) 



where the constant of proportionality depends on ro . Indeed, first note that for 
any < A < ( r °2 d d ' > ^ , it is easy to see that 

Q d n {\x"\ = X} = ( -^-x I x ew d 

which has Hausdorff dimension n — 1 , with surface measure 

n-l 



yn-l 



Then by the coarea formula, 



Q d \c = I U n ~ x {Qd n {\x"\ = \})d\ 
Jo 



\ d J 

-k 

d n-k-l 



X n ~ dX 



o 



as claimed. 
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Step 3 

We now claim that for each d > sufficiently small, we have 

x-x<u(x) (12) 

for x G dQd and x G Qd- Since u > and the function y t— > x ■ y is convex, by 
the definition of Qd we can see that it is sufficient to show (|T2"j) for all x G VV^ 
(in the case of a more general cost function, the convexity of y t— > x ■ y will be 
replaced by Loeper's maximum principle, Lemma 12.81) . 
Fix x e W d - As 

dQ d C W d , Pd U{i£ E™ | \x'\ = /(d)} U dK Pd , 

there are three cases for x G dQd- 

fie W d , w naQ d ; 

2. x £ aQ d and |.t'| = /(d) ; 

3. x g a/c Pd n aQ rf . 

Cose 1: First, suppose a; G Wd Pd C]dQd- Then since tq < 1 and by the definition 
of /(d), 

X ■ X = X ■ X + X ■ X 

< \x'\\x'\ + \x"\\x"\ 

< r /(d) 

< u(x). 

Case 2: Next, suppose that x G dQd with \x'\ = /(d). Notice that by assump- 
tion, G du(0). Thus, we have 

x-x< \x'\\x'\ + \x"\\x"\ 

< r f(d) 
= r \x'\ 

rox' 
- X ' \x'\ 

< u{x). 

Case 3: Finally, suppose that x G dK, Pd CidQd- Since the angles pd and are 
complementary, by the definitions of JC Pd and JCd we have x ■ x" = x" ■ x" < 0. 
Therefore, 

x • a; = x ■ x' + x ■ x" 

< x ■ x' 

< u{x). (13) 
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In the last line, we used the fact that x' £ du(0) for \x'\ < r . With this, we 
have shown the inequality (fT2]) for all x £ dQd and x £ Wd, thus also for all 
x £ Q d - 

Step 4 

We now claim that Qd C du(Qd), which follows from the convexity of u by 
showing that the function y i— > u(y) — y ■ x has a local minimum inside Qd for 
each x £ Qd- If this does not hold, then, it implies x ■ x < u(x) for all x £ Qd- 
However, by the calculations leading to (|13|) we see that x ■ x < u(x) for x 
in some small neighborhood of outside the cone fC Pd thus outside Qd- Since 
u(0) = 0, these together shows that is a local minimum, verifying the claim. 

Step 5 

Finally, the inclusion Qd C du(Qd) from the previous paragraph will lead to 
a contradiction. First note that since £ (sptj/) lnt , this implies if d > 
is sufficiently small we have Qd C (sptv) mt . Then from ([5]) and considering 
the volumes of Qd and Qd, this implies for all d > sufficiently small (also 
recall (HI])) 

d^-(^f)~ k \Q d \c 

< \du(Q d )nsptv\c 

< \Q d \c 

< f{d) k d n ~ k 1 

and after rearranging we obtain 

However, since n — 2k < by assumption and f(d)/d is bounded from 
above, we see that (f(d)/d) n ~ 2k is bounded below away from zero, which is a 
contradiction as d — > 0. □ 

4 General costs with MTW condition — set up 
for the proof of Theorem 11.11 

In this section, we prove some preliminary results in anticipation of the analogue 
of Step 1 (from the proof of Theorem l3.ip . for when we undertake the proof of 
Theorem O 

When the Euclidean cost function is replaced by a more general cost function 
c satisfying (|A0|) — (|A2|) and (|MTW[) . we may adapt the geometric argument used 
in the previous section to prove a result similar to Theorem 13.11 However, the 
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nonlinearity introduced by the cost function can cause a number of obstructions, 
requiring additional technical details. 

First, recall that as a preliminary step in the proof of Theorem 13.11 above, 
we "shifted" from the initial point Xo where the Theorem was assumed to be 
violated, to a point in the contact set {u = 0} which possessed the exterior 
sphere property In the Euclidean cost case we almost immediately obtained 
that the subdiffcrential of u at the original point was also contained in the 
subdifferential at the "shifted point," however for a general cost function we 
must carefully use Lemma 12751 to obtain such a containment. This is undertaken 
in Lemma |4. II below. 

Lemma 4.1. Suppose that Q and Cl are c-convex with respect to each other, 
c satisfies 4jM \i -(|A"2 ]l and (|MTW|) . and u is a c-convex function on Q. Let 
X{) G f2, and xq G d c u(xo), and define the contact set 

So := {x G fi | u(x) = —c(x, xq) + c(x , x~o) + u(x )}- 

Also suppose —Dc(xo,Xo) is contained in the relative interior of the convex set 
[d c u(xo)] X() - Then, for any point x e in So, we have 

d c u(x ) C d c u(x e ). 

Proof. Fix a point x e G So and some x G d c u{xo)- We first fix notation for the 
representations of xo, x e and xo in the cotangent spaces at xo and xo by writing 

po : = -Dc(x ,x ) S Mx , 
p e : = -Dc(x e ,x ) G [Q] So , 
po : = -Dc(x 0l Xo) G [n] xo . 

With this notation and using the definitions in © and © for c and u, we 
easily see that 

u > u{po), 
[So] XQ [fi] So I fi(p) = 5(Po)}, 

and the c-function — c(-, a;) + c(po, + u(po) is supporting to u from below at 
Po- We also define the sublevel set 

■Hi ■= {P G [«] So I -£(P, + c{p , x) < 0} 
and the (affine) half-space 

n x : = {p G T| fi | (p - po, -£>c(po, *)) < 0}. 

Now, by Lemma 12.111 (and since u > u(po) everywhere by construction), the 
sets [So] XQ and H x are convex subsets of T| Q 0. Thus by (|A1[) and (|A2|) . and by 
a differentiation of c at po it can be seen that T-L x is a supporting halfspace to 
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the set H% at p (i.e. U% C H s C T| Q Q and p e n dU s ). Additionally, 
since p e G [So] 2 we see that 

u(po) = u(Pe) > ~c{Pe, x) + 5(p , x) + u(po), 

and in particular 

Pe e vr;. 

Now, since po is in the relative interior of [d c u(xo)] x by assumption, for some 
S > we can find a c-segment x(t), defined for — 5 < t < 1 with respect to xq 
with a; at one end and xq in the interior of the curve, namely, there exists a line 
segment t G [—6, 1] t— > p(i) G [c^m^o)]^ such that 

x(t) := exp z pQ {p{t)) G d c u{x„), Vi G [-5,1]; 
xo = a;(0) and x = x(l). 

Then by the same argument as above, 

u(po) > -c(p e , x(-6)) + c(p , 5(-<5)) + u(p ) 

hence 

p e G H c x[ _ s) := {p G [fi] iQ I -c(p,5(-5)) + c(p ,5(-5)) < 0}. 
However, again by Lemma 12.111 the affine halfspace 

Hx(-s) : = {pe T? £l | (p - p Q , -Dc(p , x(-S))) < 0} 

= {pGT| o n| (p-j>o, [-CDcCso.xo)]- 1 ^-^-^))) <0} 

= {p G T| Q | (p - po, -Dcfo, sc)) > 0}. 

is a nondegenerate supporting halfspace for at po, which implies that in 

particular we must have 

Pe G {p G T| o f2 I (p - po, -Dc{p , x)) > 0}. 

Since this is the opposite halfspace to which is the supporting halfspace to 
HZ, we must actually have 

~c(p e ,x) +c(p a ,x) = 

hence 

-c(p e , x) + c(p , X) + u(p e ) = t2(p e ), 

-c(p, 5) + c(p , x) + u(p e ) = -c(p, x) + c(p , x) + u(p ) 

< u(p), Vp G [0] 5o 

which implies 

x G d c u{x e ). 

Since a; G d c u(xo) was arbitrary, this finishes the proof. □ 
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Now, we will find such a "shifted point" discussed in the preceding- 

paragraph of Lemma 14.11 Recall that in the proof of the Euclidean case (see 
Thcorcm l3.1[) . we relied heavily on the various relationships between the (n—k)- 
dimcnsional affinc set ^ containing the contact set, and the subdiffcrcntial 
du(x e ). However, in the general case the contact set between u and a supporting 
c-function, and the c-subdiffcrcntial of u at a point are contained in different 
domains. Thus we must define two systems of coordinates (actually coordinates 
on the corresponding cotangent spaces), in such a way that an argument similar 
to that of the proof of Theorem 13.11 can be applied. We step through this 
construction in Lemma l4.2l below. The statement is rather long since we include 
all the relevant assumptions. 

Lemma 4.2. Suppose that CI C M, Cl C M , c, fi, and v satisfy the same 
conditions as in Theorem \l.l\ 

Now suppose that xo G Cl is a point such that d c u(xo) H (spt v) mt ^ with 

k := affdim ([d c u{x )} Xo ) = max{affdim ([d c u(x)] x ) \ d c u(x) n (spti/)™* ^ 0}, 

and assume k > 1. Choose a point xo G Cl such that — Dc(xq,xq) G [fj] , is 
contained in the relative interior of the convex set [d c u(xo)] Xo . Finally, define 
the contact set 

S ■ = {x G ft | u(x) = —c(x, x ) + c(x ,xo) + u(x )}- 

Then, there exists a point x e G So, a basis {ei}™ =1 onT* Q, and an orthonormal 
basis {ei}™ =1 on Tj£ Q, with coordinate systems on T* Q and T£ Q centered at 
the points 

po ■= —Dc(x e ,xo) G T* CI and p e := —Dc(x e ,x ) G T^Cl (14) 

with coordinate directions given by {ej}™ =1 and {ei}™ =1 , respectively, and positive 
numbers Ro and ro with the following properties (Recall Remark \2.1\ on notation 
using coordinates): 

[So] So C{pe [Cl]- Xo c TICl \p' = 0, \p" + Roe n \ < R a }, (15) 

P-P = (p - Pe, [-DDc(x e ,x )}~ 1 (p - po)), 

VpGT^n, peT* e Ct, (16) 
fl*(0) : = {peT*n | \p'\ <r , p" = 0} 

C [d c u(x e )] Xe , (17) 
Hl ro - ei c{ P e [Cl] xo | ± P l < 0}, VI < i < k. (18) 

Here, Hp is defined by (|20p below. 

Proof. Since CI is strongly c-convex with respect to spt v, and 9 c M(xo)n(spt z/ mt ^ 
0) by © we may use the boundary- not-to- interior lemma in [131 Theorem 5.1 
(b)] to see that the closure of So is contained in the open set CI, i.e. Sq 1 C CI. 
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This shows in fact, So = Sq. Thus, [Sq] s is closed and also by Lemma B.lll it 
is convex. Therefore, there exists a point p e S d [So} £ with the exterior sphere 

condition and withp e S P]" 1 ', which will lead to (|15l) after choosing appropriate 
coordinates below. Define 

x e := exp% Q (p e ) e dS . 

Now, since 9 c w,(a;o) (~l (sptz/) lnt ^ by assumption and d c u(xo) C d c u{x e ) by 
Lemma 14.11 above, we find d c u(x e ) fl (spt ^) lnt ^ as well. The inclusion 
d c u(xo) C 9 c u(x e ), and the maximality of k implies 

affdim ([d c u(x e )] x ) = k. 

Also, since the map —Dc{x e ,exp% {-)) : [d c u(xo)] Xo — > [d c u{x e )] x is a diffeo- 
morphism, the point 

Po ■= -Dc(x e ,x ) 

is contained in the relative interior of the convex set [d c u(x e )] x . 

We must now determine an (n — fc)-dimensional affine space containing the 
contact set, as in ((6]) (we must, however, point out that the (n — fc)-dimensional 
set determined here is only for the purpose of defining coordinates, in the proof 
of the actual theorem we will be forced to choose a curved analogue of (O ) ■ To 
this end, define for each p S [d c u(x e )] x the point 

x p : = exp c Xe {p) 6 d c u(x e ). 

We use the definitions ^ and ([3]) for c and u (modified with respect to the 
fixed point xq), and consider the c- functions associated to Xp, 

rhp{p) : = -c(p, xp) + c(p e , %) + u(p e ), (19) 

the c-sublevel sets 

Hp : = {pe [Q] So | m p {p) < u(p e )}, (20) 

and affine half-spaces 

Hp : = { P e T! o n | (p - p e , -Dm(p„)> < 0}. 

As in the proof of Lemma 14.11 above, by Lemma 12.111 each set Hp is convex in 
T| o fi, and by (|A1[) and (|A2|) . Hp is a supporting halfspace for the set Hp at p e 
for each p£ [f2] . Since 

[S„W = {pe I 2(p) = fi(Pe)} 
and each is supporting to u, we find that 

pe[a c u(x c )] Xc 

c n ^- 
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Since po is in the relative interior of the fc-dimensional set [d c u(x e )] x , it follows 
that A x ~ k is an (n — fc)-dimensional affine set. 

We now choose an orthonormal basis (in the Riemannian metric g So on Tg f2) 
{ e i}™=i 01 T^ Cl such that {ei}™ =fe+1 span the linear subspace (A x ~ k — p e ) ■ Also, 
since p e satisfies the exterior sphere property in [So] x , we may choose the basis 
and coordinate directions so that (|15p is satisfied for some Rq > 0. 

Next, we must determine the appropriate sense of "orthogonality" between 
the sets A x ~ k defined above, and [d c u(x e )] x . Choose a collection of vectors 
in T* given by the dual relations 

Sij = (e t , [-DDc{x e ,x )\~ X ej) 

for each 1 < i, j < n. Since {e.i}™ =1 is actually a basis by (| A2|) . we may 
also define a coordinate system (denoted by p) centered at po with positive 
coordinate directions given by {ei}" =1 . Then, for any p e T^Q and p £ ft, 
we calculate 

n 

(p - p ei [-DDc(x e ,x )]~ 1 (p - p a )) = ^2 (p %e h[-DDc(x e ,x Q )]~ 1 p :i ej) 

n 

= ^2 P^ieui-DDcixe,^)]- 1 ^) 

i,j=l 
= P-P 

and hence obtain (fT6)l . We easily see from the definition of A™~ k that {ei} k =1 
spans the linear subspace (aff ([d c u(x e )] x ) — po)- Since p is in the relative 
interior of the convex set [d c u(x e )] x , this implies that we obtain (jTT)) for some 
r > 0. Finally, using (fTl)|) we calculate for any fixed 1 < i < k that 

%±r e< C ^roet HI [fi] io 

n 

= {pe [fi] So I ±n,^^'(e;, [-DDcfaxo)]- 1 ^) < 0} 

3=1 

= { P e [fi] Xo | < 0} 

to obtain the final property (fT5|) . □ 

5 Proof of Theorem 11.11 

We are finally ready to give the proof of Theorem 11.11 The results in Sec- 
tion 2] allow us to make geometric constructions similar to those in the proof of 
Theorem 13. II Still, the nonlinearity of the cost function produces a number of 
additional difficulties. 

For one, the set A x ~ k from the proof of Lemma [4721 above is not sufficiently 
geometrically motivated enough to fill the role of A n ~ k from the proof of Theo- 
rem [3J] Instead, we must use A™T k (defined in (|22|) below) which is no longer 
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an (n — /s)-dimensional affinc set, but a locally smooth (n — fc)-dimensional sub- 
manifold of O. We must carefully consider its relation to the coordinate system 
we defined above in Lemma 14.21 Also, with a general cost c, we cannot expect 
the vanishing of terms due to orthogonality that we enjoyed at various points 
in the proof of Theorem 13. II Instead, we must carefully keep track of the size of 
such error terms and show that they can be ultimately controlled. A side effect 
of this last difficulty is that we must alter the definition of f(d) to account for 
the rate of decay of u(x) as x approaches the contact set. 

Remark 5.1. In the interest of readability, we make the following comment 
concerning constants in the following proof. We will absorb all constants that 
depend only on the cost function, the domains f2 and Cl, and the constants n, 
fc, ro, i?o which arc fixed at the beginning of the proof, into the single constant 
AI. There are three constants j3, Co, and C\ which are introduced and left "to 
be determined" , we keep careful track of these constants throughout the proof. 
These constants are necessary in order to control various error terms that arise, 
due to the nonlinearity of the cost function. Eventually we will determine that it 
is sufficient to take /3, Co, and C\ to be small, depending on M but independent 
of the parameter d. 

5.1 The Proof 

The proof is outlined similarly as the one for Theorem 13. II 
Proof of Theorem 

Step 

Suppose that the theorem does not hold. Since d c u(x) = exp%(du(x)) for each 
x G by Corollary 12. 101 we see that this implies 

k : = affdim ([d c u(x )} Xo ) 

= max {affdim ([d c u(x)] x ) \x eQ and d c u(x) n (spt v) int ^ 0} 

for some Xq £ f2 with k > ^. Note that by ([1} we must have k < n— 1. 
Step 1 

We now apply Lemma 14.21 above to obtain points x e € f2, xq £ Cl (and corre- 
sponding p e and po, see (fT4)) ). coordinate bases on T^ g Cl and T* e Q, denoted by 
{ej}™=i and {e;}™ =1 respectively, and positive numbers Ro and ro. Instead of 
working with the original c and u, we will use in the following, the modified 
functions c and u with respect to the point xq (see the definitions @ and ([3]) 
). By adding a constant we may assume u(p e ) = 0, and by translations of the 
domains assume that both 

p e and pa = 



19 



(hence both coordinate systems on [f2] s and [f2] are centered at the origin). 
In particular, this implies 

u(0) = 0, u > 0, and e d £ u(0). 

Also, by an abuse of notation we will identify a point p £ [Q] £ or p £ 
[n] with the coordinates given by {e;} or {ei}. In particular, note that by 
property Q16p we have 

p-p={p,[-DDc(x e ,xoT 1 p) (21) 

for any p £ [^] So an d p£ [fi] i while by (|A1|) and (|A2[) . the length of any 
vector p in the Riemannian metric on T* Q is comparable to the Euclidean 
length \p\ measured in these coordinates. In the remainder of the proof, it is 
helpful to keep in mind that p' £ dcu{0) for any \p'\ < r^. 



Step 2 

Again as in the proof of Theorem 13.11 we aim to construct two families of sets 
Qd C [0] s and Qd C [f2] depending on a small parameter d: the ultimate aim 
is to apply inequality ([T]) and obtain a contradiction in the volume comparison 
as d approaches zero. First, define the set (with notation as in Lemma 14.21 and 
(fH)j) above) 

k 

At k :=C\{P 6[ fi UI^(p)=0} (22) 

i=l 

(compare with (|6])): by transversality of each level set {p £ | "V s;(p) = 0}, 

is a locally n — fc-dimcnsional submanifold, but not necessarily flat. 
Choose 

7T (Z 

Pd " 2 + 4i^' 

for small d > 0, and define the n-dimcnsional cone in [£l] s with angle pd by 
£ Pd ■= {p e [^] So I |p| cospd < p n }. 

Also, define 

w d , Pd := { P £[n} Sa \\p"\ = d}rMC Pd . 

Then from the same calculations leading up to ((9]) we obtain by (|15|) that 

w d , Pd n [S ] So = 0. 

In particular, 

■ti(p)>0 for peW d , Pd . 

Since f2 is bounded, as in the proof of Theorem 13. 1[ we can define the strictly 
positive quantity 



M 



d 2 , min u(p) 
peW d , Pd ' 



> 
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for each d > sufficiently small. The definition of f(d) differs from the proof 
of Theorem 13.11 the reason being that we cannot rely on orthogonality to fully 
eliminate all of the extraneous terms as we did for the Euclidean cost. 

Before we construct the families of sets Qd and Qd to be used in the volume 
comparison argument, we first give an appropriate parametrization of the pos- 
sibly curved set A™~ k . Recall that we are equating a point in [f2] 2 with a point 
in W l by its coordinates with respect to the basis {e^}. Define the mapping 
F : [fi]^ -> R k by 

F{p) := (m roSl 0), . . . , m roSk (p)), 

and ip : M$ -> t" by 

tl>(p) :=(F(p),p"). 

Note that by construction F(0) = 0, thus ip(0) = 0. By the inverse function 
theorem, ip is a diffeomorphism for p in some small neighborhood (depending 
only on c) of the origin (recall, which is contained in the interior of ), and 
we can write the inverse as 

^-\ P ) = (G(p),p") 

for some map G : R" -> {p E [Cl]- o \ p" = 0}. Notice that F^^ip")) = thus, 

^ip") e Ar k - (23) 

The vector p — ip'~ 1 (p") is in the span of the first k basis vectors ej, 1 < i < k, 
and gives the displacement of p from the set A™~ k . This vector will be used 
crucially in Step 3. The size of this displacement is 

\p-^- 1 {p")\ = \p' -G{p")\ 

= \G^(p))~G{p")\ 
~m \tP(p)-p"\ 
= \F( P )\ 

i.e., 

M- l \F(p)\ < \p-^- l (p")\ < M\F(p)\. (24) 

Also note that since F and G have Lipschitz constants depending only on c and 
spt v, we can calculate 

\p\ < \p'\ + \p"\ 
= \G{i>(p))\ + \p"\ 

< \G{i>(p)) - G(p")\ + \G( P ") - G(0)\ + \ P "\ 
<M(\iip)-p"\ + \p"\) + \p"\ 

< M(\F(p)\ + \p"\) (25) 
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while since "0(0) = we find 

N» -1 (P)I ~M bl- (26) 
Now, define the desired family of sets in f2 by 

Q d := { P e M So I \F(p)\ < /(d), |p"| < d} n AC W . 
If d is sufficiently small, the set 

{peft I |f(p)| < /(d), <d} 

is diffcomorphic to the set 

{x g M' 1 | < M/(d), < d}, 

hence we have the volume bound 

\Qd\c<f{d) k d n - k . (27) 

To define the other family of sets Qd in 0, consider 

%p d -./ 2 ) [n] a< | |p"|cos(/3( Pd -7r/2)) <p»}, 

cylinder whose base is an (n— fc)-dimensional cone in [fi] ^ with angle /3(/9<i — 7r/2). 
Using this we define the set (compare with (|ll)p) 

W d := iCp {pd _ w/2) n{pe | P' G conv{0, Code, | 1 < i < fc}, = Ci^} 

(28) 

for appropriate small positive constants /3, Co, and Ci to be determined, inde- 
pendent of d. We now define the target shape Qd as the cone over Wd with 
vertex 0, that is, 

j>ew d 

where [0,p] denotes the line segment (in the convex set [Q] ) from to p. By 
a calculation similar to II H we see Qd has the volume bound 



n— k 

-k-lik 



d™-*- 1 ^ ~\Q d \c (29) 



Step 3 



This step will be the most involved. We now claim that for some appropriate 
choice of f3, Co and C\ above (independent of d), for each d > sufficiently 
small we have 

l fnp(j ) ) < u(p) for every p E dQd as long as p € Qd (30) 
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(compare with (fT2"j) ). Notice that u > and mo = (because c(-,Xo) = 0), 
moreover, rhp satisfies Loeper's maximum principle (Lemma I2.8|) along each 
line segment in [0] x . Therefore, by the definition of Qd we sec that it is 

sufficient to show ([3TJ]) for each p £ Wd- 
Fix p g W d - As 

9Q d c W d , Pd u{ P e [fi] So | = /(d)} u 9/C Pd , 

there are three cases for p S <9Qd: 

1. pe W d , Pd n 9Q d ; 

2. pedQ d and |F(p)| = /(d); 

3. p 6 d/C Pd n aQ d . 

Case J: First suppose that p £ Wd, Pd n we then calculate for sufficiently 
small d, 

fhpip) = [fhpip) - fhp(p)] + [rhp>{p) - rh pl (^ 1 (p // ))} + rhp,{^ l {p")) 
= I + II + m p ,{i>- 1 (p")). 

Then, using (|25|) . we find, 

I = -c(p, %) + c(0, x p ) - (-c(p, %) + c(0, %)) 
<M|p||p-p'| 

<M(|F(p)| + |j/W'l 

< M(f(d) + d) ( since p G Vfc, |p"| < ) 

< Md/(d) 

for some M > depending only on bounds on derivatives of c. Similarly, we 
see by using that 

7/ = -c(p,%<) +c(-0~ 1 (p"),%) 

< M|p- V _1 (p")IIp'I ( sm ce c(-,x ) = 0) 
<M\F(p)\\p'\ 

< MCodf(d). 

Finally, if d is small enough, then for p G Wd, p' is in the convex hull conv{roe"i 
1 < i < k}. 

< max {"V e i (V , ~ 1 (p")} 

1<?'</c 

= 
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by Lemma |2~51 and (|2"5]l . Hence all together for small constants Co, C\ (depend- 
ing only on M) we obtain 

m p {p)<M[C a d+C 1 ]f{d) 

<f{d) 
< u(p) 

where the last line follows from the definition of f(d), since p £ Wd, Pd - The 
claim (|30)) is verified in this first case. 

Case 2: Next, suppose that p £ dQd, and \F(p)\ = f{d). This case is subtle. In 
the following, the vector p — will play an important role. Recall that 

it is contained in the span of the first k basis vectors, ej, 1 < i < k. Thus by 
the relations (fT6|) and using (|A2[) . for some choice of either plus or minus, and 
some index 1 < jo < fc, we find that 

±(p-rV))-e i0 > j^b-rVOi. (si) 

where M > depends only on c and k. At the same time, we claim that 

m p , (p) < fh± roSjo (p), V P 6 Xr fe - (32) 

Indeed, note that p' £ conv {0, r ei 1 < i < k} if d is small enough. Thus by 
applying Locper's maximum principle (Lemma 12. 8|) multiple times we see that 

rhpip) < max{m (p), m ro g io (p)} 
= 

= rn ro g Jo (p) 

if p £ (using that mo = 0). On the other hand, by property (jT5)) we have 

{p £ [0]- o | m_ roe - 3 . < 0} c {p £ [fi] 2o I V < 0} 
hence by the continuity of c, 

AT k c {p £ [n] So | p jo < o} £ {p £ [n] So | < m_ roe - jo }, 

and we obtain the claimed inequality (|32[) . It can be seen that we may assume 
the choice of vector is +e~j without any loss of further generality, so we shall 
do so for the remainder of this case. 

Now our goal in this case is to show the inequality fhp{p) < 7h rogj (p). 
Indeed, roe~j a £ <9gu(0) implies rh ragjo < u, which then easily leads to the desired 
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inequality ([50)) . We assume in the following d is sufficiently small. First write 

fhpijp) = [rh p (p) - fhp>{p)] 

+ [rh p ,(p) -rh pl {^- l {p"))] 

+ [-m ro s JO (p) + < e - {^{p"))] 

+ [m pl ^- 1 (p"))-rh roS J^- 1 (p"))} 

+ ™r e jo (P) 

= 1 + 11 + III + IV + rh roSjo (p) 
Here, similarly to Case 1 , we see 

/ < Md/(d). 

To deal with the term 77, we first make an auxiliary calculation, which will 
be of use in a number of places later. Fix any pi, p 2 G [^ls an d P & [&] x ■ 
Then using a Taylor expansion argument along with (| A2|) , ((4]) , and (fT6|) results 
in the bound 

-c(p 2 ,Xp) +c(pi,%) 

< (-£>c(pi,%)),p 2 - Pi) + A/|p||p 2 - Pi I 2 

< (-7>c(0,%)),p 2 - pi) + M\p\\ Pl \\ P2 - Pl \ +M\p\\p2 - pi | 2 

= p-(p 2 -Pi) + MIpIMIps -pil +M\p\\p 2 - Pl \ 2 . (33) 

By applying (f33"f to the term 7/ and using the estimates (f2"5j) and (|2"6]) . we 
obtain 

77 = -c(p,%>) + c(0,%) - (-c(V'~ 1 (p"),%) + c(0,%)) 

< p' ■ (p - rV')) + m^-^jiip - ^-V)l + Mb - ^ 1 (p")I 2 

< \p'\\F(p)\ + M\p\\F(p)\ + M\F(p)\ 2 

< C Mdf(d) + M(f(d) + d)f(d) + Mf{d) 2 

(using pe Qd,P& Qd) 

To deal with the term 17/, first we note that since fho = 0, by (|A2|) we will 
have 

\Dm p \ < M\p\. 
Hence by {33]), ||25J), and {25) again, 

^ = -CO/*" 1 (p"), ^r e J0 ) + C(p, X rgSj0 ) 

< -r e j0 ■ ( P ~ ^ V')) + M\i>-\p")\\p - i>-\p")\ + M\p - ^ V')| 2 

< — M _1 |p - V>- V')l + M\p\\p - i)-\p")\ + M\p - t/T V)| 2 

(the first term by [[31])) 

< -M-V(d) + M(/(d) + d)/(d) + A//(d) 2 

(using pe Qd, P ^ Qd and the condition |F(p)| = f(d)). 
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The key here is that the first term —M~ 1 f(d) in the last line above decays 
slower than all other terms in I, II and 77/, if C\ is chosen to be small. In 
particular, for sufficiently small fixed C\ and for small d, we see 

I + 11 + III < 0. 

Finally, since ip" 1 (p") G A™~ k wc have IV < by (J32J). All together, this shows 
that rhp{p) < m roSjg (p) as desired, and we obtain (f5U|) in this case. 

Case 3: 

Finally, suppose that p € dIC Pd H dQd- This case is subtle as well. First, 
notice that for small universal constants /3 and Co, and for small d > 0, for 
p G fC Pd and p G Wd, the angle between p and p" = {0,p") is greater than 
7r/2 + "fd with 

Note that by applying similar calculations as (|33|) but in the variables p, we 
obtain 

[-c(p, x P2 ) + c(0, % 2 )] - [-c(p, x Pl ) + c(0, % )] 

<p-(p 2 - pi) + M|p||pi||p 2 ~Pi\ + M\p\\p 2 -pi\ 2 
Combining this with the lower bound on the angle above, we obtain 

= [-c(p, x p ) + c(0, %)] - [-c(p, x p > ) + c(0, xp> )] 
<p-p" + M\p\\p'\\p"\ + M\p\\p"\ 2 

< -\p\\p"\sm(jd) + M\p\\p"\\p\ 

(from the lower bound on the angle between p and p") 
<-M(l-p)\p\\p"\d + M\p\\p"\\p\ 
= \ P \\p"\(-M(l-f3)d + M\p\) 
<(j3 + C + C 1 -l)M\p\\p"\d 

(since peW d and |p| < C d + df(d)/d < (C + C^d ) 

< (34) 

for sufficiently small constants /3, Co, C\ (independent of d), and for sufficiently 
small d. As a result, 

fhpip) < mp'(p) < u(p) 

where we have used the fact that p' <G 9gu(0) for \p'\ < r$. This shows (|30[) for 
the final case. 

With this we have shown the inequality (|30|) for all p G dQd and p G VVd, 
since it clearly holds when p = the inequality holds for all p G Qd by 
Lemma 12.81 
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Step 4 



We now claim that Qd C d c u(Qd) for d sufficiently small. If we can first show 
that u — rhp has a local minimum in the set Qd for each p £ Qd, we may then 
apply Corollary 12.101 to obtain the desired inclusion. Indeed, if there exists 
a point in the interior of Qd where u < rhp we immediately obtain a local 
minimum from ([30]) . If this does not hold, we see by the calculations leading 
to (|34p that rhp < u on a small neighborhood of 0. Since u(0) = = rhp(0) by 
construction, we see that is a local minimum of u — rhp. 



Step 5 

Finally, by using (fTJ) in conjunction with the volume bounds ([27]) and ([29]). we 
see (for all d > sufficiently small) 



7i— k 

i—k—1 ik 



< \d c u(Qd)\c 

< \Qd\c 

< f(d) k d n ~ k . 

However, since n — 2k < by assumption and f(d)/d < d by definition, we 
see that (/ '(d) j d) n ~ 2k is bounded below away from zero, so after rearranging we 
obtain 

l<(^)"~ 2fe <^. 
which is a contradiction as d — > 0. This completes the proof of Theorem 1 1. 1 1 □ 



Appendices 

A An example of a purely (IMTWf ) cost 

In this section, we give an example of a cost function which satisfies (|MTW[) , but 
not ( MT W+ I and (|NNCC[) . This example originally appears in [32] where it is 



stated to satisfy (]MTW[) but not QMTW+fl . It was communicated to the authors 



by Neil Trudinger [30] that the cost function actually fails to satisfy ([NNCCI) 
as well. The proof given in this appendix is only for the reader's convenience, 
since the detailed proof docs not seem to appear in the literature. 

Lemma A.l. Define 

c(x, x) :~ \x — x\~ 2 

for x, x e R n . Then, c satisfies conditions (fX0]) - ([X2]) and ([MTWP . but does 
not satisfy ( |MTW + [ ) and (jNNCC[) on any subdomain o/M™ x E" whose closure 
does not contain the diagonal {(x,x) | x = x}. 
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Proof. It is clear that c satisfies (|A0[) on any such a subdomain of 1" x R n . 
Assume throughout that x ^ x. We can calculate, 



-Dc(x, x) = 2{x — x)\x — x\~ 



it is clear that for a fixed x, this is an injective mapping in x, and by the 
symmetry of the two variables, we see that (|A1I) is satisfied. Additionally, 



DiDjc(x,x) = —D^c(x,x) 



hence 



2\x — x\ 4 [Sij — 4(xi — Xi)(xj — Xj)\x — x\ 2 
4(x — x) ® (x — x) 



(35) 



dct DDc{x, x) = 2 n \x- 5|- 4,l dct Id 



= -Z-2 n \x~x\- An ^0, 

hence we obtain (|A2|) . 

We will now verify condition (|MTW|) by a lengthy, but routine calclation. 
Note that an alternate way of writing condition (|MTW[) is by the inequality 

for x € R™, p 7^ 0, and (V, r/) = 0, where 

A i:j (x,p) := -D 2 i:j c(x,expl{p)). 

Now for p £ R", since by definition —Dc(x, exp^.(p)) — p, by taking the absolute 
value of both sides of this expression we find 

2-*\p\* = \x- exp c x {p)\-\ 

which leads to the formula 

x-expl{p) = 2*\p\-*p. 
Plugging this expression into ([3"5j) , we find 

Aijfap) = 2~*bl*(*« ~^Pj\p\~ 2 )- 
Differentiating once in p, 

= 2-i[-pk]p\~*(6 ij - ipiPj\p\~ 2 ) - i\p\Hs ikPj \p\- 2 + Pi5 jk \p\- 2 - 2p i p jPk \p\~ i )} 



2'W 



1 2 

^hjPk - SikPj -PiSjk + ^PiPjPk\p\ 



28 



and differentiating in p again, 



D l k P, A ij( x >P) 



2 3 ( - 3 )Pl\p[ 



1 2 

-Sijpk - S ik pj -PiSjk + -^PiPjPk\p\~ 



2*\p\~ 



-^SijSki — SikSji — SuSjk 



2 2 2 4 



2§ i 



-g%j»fcP/bl 2 + 25 ik PjPi\p\ 2 + 2p l S jk pi\p\ 2 



+25uPjp k \p\~ 2 + 2piSjipk\p\~ 2 + 2p t pj5ki\p\~ 2 - ^PiPjPkP^pf 1 } 
Hence, for any V, r\ £ 1" we find 
DlA^p^V^M 



2§ 



■w 



l^fl'^ 2 + 2{v > v)i W\ ,v){ V\ ,v) + 2{v,v){ W\ ,v){ ¥\ ,ll} 



-l(^,v) 2 (^,v) 2 + \v\ 2 \v\ 2 -mv) 2 -Wv) 2 

+2(F, ?? )(^, ?7 )(^,^) + 2^, ?? )(^,F}(^, ?7 )+2| ?7 | 2 (^,F} 2 -4(^,F} 2 ( 



4w 



|y|^ + 2b| 2 <Z,F> 2 - = |V| a <iW - ^(f-.y^Jl 



2/ P 



16, p m2/ P 



-6(y ?? ) 2 + 8(y,r 7 )(^,y)(^,r ? ) 



'bl" '>l 



bl 



3 bl 



'bl 



We are now ready to verify (|MTW[) . Let us write 

<1 



P_ 

bl' 



and assume that \V\ = \r}\ = 1, with (V^, 77) = 0. By solving the constrained 
optimization problem of maximizing the function q 1— > (q, r\) 2 + (q, V) 2 subject 
to |<?| 2 = 1 in E™, we see that 



\(q,v)\ < 
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and thus for (V, 77) = 0, 



2~ 
"3 

2l 



> T w 



1 + 29 2 



(1 



16 



2 (l-0 2 ) 



2- 



i6<r 



This expression attains the minimum value of zero for 9 G [0,1] at = 1/2, 
hence we see that c satisfies (jMTWI) . 

Next, by the calculation above, we see that if we choose p £ M. n \ {0} and 

(V,r?) = with |V| = \r]\ = 1 so that 9 = ±± and (g,r/) = Jl-\ = # (for 



example, V = ei, ry = e2 and p 



-^62), then we actually have 



hence c does not satisfy ( |MTW + 



Finally, we show that c does not satisfy (jNNCCp . Consider V = A77 with 
I?/ 1 = 1 and any A 7^ 0. Then for any p that is a nonzero multiple of rj, 

Dl kV A l2 {x,p)V'V^ m 



21 1 1-^2 
= T bl » A 2 



2" 



H^A 2 



< 



finishing the proof. 



□ 



References 

[1] Yann Brenier. Polar factorization and monotone rearrangement of vector- 
valued functions. Comm. Pure Appl. Math., 44(4):375-417, 1991. 

[2] Luis A. Caffarclli. A localization property of viscosity solutions to the 
Monge- Ampere equation and their strict convexity. Ann. of Math. (2), 
131(1):129-134, 1990. 

[3] Luis A. Caffarclli. Boundary regularity of maps with convex potentials. 
Comm. Pure Appl. Math., 45(9):1141-1151, 1992. 

[4] Luis A. Caffarelli. The regularity of mappings with a convex potential. J. 
Amer. Math. Soc., 5(1):99-104, 1992. 



30 



[5] Luis A. Caffarelli. A note on the degeneracy of convex solutions to Monge 
Ampere equation. Comm. Partial Differential Equations, 18(7-8):1213- 
1217, 1993. 

[6] Luis A. Caffarelli. Allocation maps with general cost functions. In Partial 
differential equations and applications, volume 177 of Lecture Notes in Pure 
and Appl. Math., pages 29-35. Dekker, New York, 1996. 

[7] Piermarco Cannarsa and Yifeng Yu. Singular dynamics for semiconcave 
functions. J. Eur. Math. Soc. (JEMS), 11(5):999-1024, 2009. 

[8] Philippe Delanoe. Classical solvability in dimension two of the sec- 
ond boundary-value problem associated with the Monge- Ampere operator. 
Ann. Inst. H. Poincare Anal. Non Lineaire, 8(5):443-457,1991. 

[9] Philippe Delanoe and Yuxin Gc. Regularity of optimal transport on com- 
pact, locally nearly spherical, manifolds. J. Peine Angew. Math., 646:65- 
115, 2010. 

[10] Philippe Delanoe and Yuxin Ge. Locally nearly spherical surfaces are 
almost-positivcly c-curved. Methods Appl. Anal., 18(3):269-302, 2011. 

[11] Alessio Figalli, Young- Heon Kim, and Robert J. McCann. Continuity 
and injectivity of optimal maps for non-negatively cross-curved costs. 
arXiv.09ll.3952, November 2009. 

[12] Alessio Figalli, Young-Hcon Kim, and Robert J. McCann. Regularity of 
optimal transport maps on multiple products of spheres. arXiv: 1006. 1957, 
June 2010. To appear in J. Eur. Math. Soc. (JEMS). 

[13] Alessio Figalli, Young-Hcon Kim, and Robert J. McCann. Holder continuity 
and injectivity of optimal maps. arXiv:1107.1014, July 2011. 

[14] Alessio Figalli, Ludovic Rifford, and Cedric Villani. Nearly round spheres 
look convex. Amer. J. Math., 134(1):109-139, 2012. 

[15] Alessio Figalli and Ludovic Rifford. Continuity of optimal transport maps 
and convexity of injectivity domains on small deformations of E> 2 . Comm. 
Pure Appl. Math., 62(12):1670-1706, 2009. 

[16] Alessio Figalli, Ludovic Rifford, and Cedric Villani. On the Ma-Trudinger- 
Wang curvature on surfaces. Gale. Var. Partial Differential Equations, 
39(3-4) :307-332, 2010. 

[17] Alessio Figalli, Ludovic Rifford, and Cedric Villani. Necessary and suffi- 
cient conditions for continuity of optimal transport maps on Riemannian 
manifolds. Tohoku Math. J. (2), 63(4):855-876, 2011. 

[18] Wilfrid Gangbo and Robert J. McCann. The geometry of optimal trans- 
portation. Acta Math., 177(2):113-161, 1996. 



31 



[19] Young-Heon Kim and Robert J. McCann. Continuity, curvature, and the 
general covariance of optimal transportation. J. Eur. Math. Soc. (JEMS), 
12(4):1009-1040, 2010. 

[20] Young- Hcon Kim and Robert J. McCann. Towards the smoothness of 
optimal maps on Riemannian submersions and Ricmannian products (of 
round spheres in particular). J. Reine Angew. Math., 664:1-27, 2012. 

[21] Jun Kitagawa. An iterative scheme for solving the optimal transportation 
problem. arXiv:1208.5172, August 2012. 

[22] Jiakun Liu. Holder regularity of optimal mappings in optimal transporta- 
tion. Gale. Var. Partial Differential Equations, 34(4):435-451, 2009. 

[23] Jiakun Liu and Neil S. Trudinger. On Pogorelov estimates for Monge- 
Ampere type equations. Discrete Contin. Dyn. Syst., 28(3):1121-1135, 
2010. 

[24] Jiakun Liu, Neil S. Trudinger, and Xu-Jia Wang. Interior C 2 a regularity for 
potential functions in optimal transportation. Comm. Partial Differential 
Equations, 35(1):165-184, 2010. 

[25] Gregoire Loeper. On the regularity of solutions of optimal transportation 
problems. Acta Math., 202(2):241-283, 2009. 

[26] Gregoire Loeper. Regularity of optimal maps on the sphere: the quadratic 
cost and the reflector antenna. Arch. Ration. Mech. Anal., 199(l):269-289, 
2011. 

[27] Gregoire Loeper and Cedric Villani. Regularity of optimal transport in 
curved geometry: the nonfocal case. Duke Math. J., 151(3):431-485, 2010. 

[28] Xi-Nan Ma, Neil S. Trudinger, and Xu-Jia Wang. Regularity of poten- 
tial functions of the optimal transportation problem. Arch. Ration. Mech. 
Anal., 177(2):151-183, 2005. 

[29] Robert J. McCann. Polar factorization of maps on Riemannian manifolds. 
Geom. Fund. Anal, ll(3):589-608, 2001. 

[30] Neil S. Trudinger. Personal communication. 

[31] Neil S. Trudinger and Xu-Jia Wang. On strict convexity and continuous 
differentiability of potential functions in optimal transportation. Arch. Ra- 
tion. Mech. Anal., 192(3):403-418, 2009. 

[32] Neil S. Trudinger and Xu-Jia Wang. On the second boundary value problem 
for Monge-Ampere type equations and optimal transportation. Ann. Sc. 
Norm. Super. Pisa CI. Set. (5), 8(1):143-174, 2009. 

[33] John Urbas. On the second boundary value problem for equations of 
Monge-Ampere type. J. Reine Angew. Math., 487:115-124, 1997. 



32 



[34] Cedric Villani. Stability of a 4th-order curvature condition arising in opti- 
mal transport theory. J. Fund. Anal., 255(9):2683-2708, 2008. 



[35] Cedric Villani. Optimal transport: Old and new, volume 338 of Grundlehren 
der Mathematischen Wissenschaften [Fundamental Principles of Mathe- 
matical Sciences]. Springer- Verlag, Berlin, 2009. 

[36] Xu-Jia Wang. Some counterexamples to the regularity of Mongc- Ampere 
equations. Proc. Amer. Math. Soc, 123(3):841-845, 1995. 



33 



